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ABSTRACT 

We compare the simultaneous 2007 space-based MOST photometry and ground-based radial velocity 
observations of the F5 star Procyon. We identify slow variations in the MOST data that are similar 
to those reported in the radial velocity (RV) time series, and confirm by comparison with the Sun that 
these variations are likely the signature of stellar activity. The MOST power spectrum yields clear 
evidence for individu al oscillation frequencies that match those found in the radial velocity data by 
iBedding et all (|2010T ). We identify the same ridges due to modes of different spherical degree in both 
datasets, but are not able to confirm a definite ridge identification using the MOST data. We measure 
the luminosity amp litude per radial mod e ^z^o.phot = 9.1 ± 0.5 ppm. Combined with the estimate 
for the RV data bv lArentoft et al.l (|2008D this gives a mean amplitude ratio of A ;=0 . P hotMi=o,RV = 
0.24 ± 0.02 ppm cm _1 s, considerably higher than exp ected from scaling relations but in reasonable 
agreement with theoretical models bv iHoucfe k (2010). We also compare the amplitude ratio as a 
function of frequency, and find that the maximum of the oscillation envelope is shifted to higher 
frequencies in photometry than in velocity. 

Subject headings: stars: individual (Procyon) — stars: oscillations — stars: activity — stars: rotation 
— techniques: photometric — techniques: radial velocities 



1. INTRODUCTION 

The detection and measurement of oscillations in stars 
provides a unique possibility to infer details about the 
physics governing their interiors. The prospect of extend- 
ing such studies from the Sun to distant stars has mo- 
tivated many observation campaigns in recent decades. 
Even with the wealth of new space-based photometry 
from CoR oT (see, e.g.. IMichel et al.l 120081) and Kepler 
(see, e.g., Gilliland et al.l I2010D there is still an impor- 
tant place for ground-based spectroscopic campaigns of 
bright nearby stars with well-known fundamental param- 
eters. Owing to its brightness (V = 0.3), proximity (d = 
3.5 pc), and membership in an astrometrically well de- 
termined binary system, the F5 sub-giant Procyon A 
(aCMi, HR2943, HD 61421) has long been considered 
a prime target for such campaigns. 

The majority of early efforts to detect oscillations in 
Procyon have relied on measuring Doppler velocities 
from a single site. The fi rst claimed detection dates 
back to IGellv et al.l (|1986f ) which, however, could not 
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be confirmed bv iLibbrechtl (|1988t ) or llnnis et all (|1991l) . 
who reported null-detections at similar sensitivity lev- 
els. With the benefit of hindsight, it now seems that 
the first detection of power excess in Procyon (and, in 
fact, any other so lar- like star than the Sun) was by 
iBrown et all (fl99lh . This was followed by numerous ob- 
serving campaigns, mostly single-site, taking advantage 
of the increasin g prec i sion of Dopple r-shift measurements 
(iMosser et all Il998t iBarban et all 119991: QVIartic et all 
1 19991 12004t iClaudi et al.l 120041: [Eggenberger et all 120041: 
iBouchv et al.H2004tlLeccia et al.ll2007f) . While alTof these 
studies revealed clear power excess in the expected fre- 
quency range of 0.5-1 mHz, a consistent determination of 
individual frequencies was hampered due to severe alias- 
ing caused by daily gaps inevitable in dual or single-site 
observations. 

The first two sets of continuous observ ations of Pro- 
cyon by the Canadian MOST satellite (| Walker et all 
12001 iMatthewsl I2007T ) in 2004 and 2005 resulted in 
null-detections, leading to the conclusion that lumi- 
nosity amplitudes in Procyon must be lower than 
15 ppm and/or the mode life t imes shorter than 2 - 
3 davs dMatthews et all l2004t iGuenther et~aT1 120071 ). 
IBedding et al.l (|2005) found these results to be com- 
patible with limits set from ground - based radial ve- 
locity observations a nd iBaudin et all (120081) confi rmed 
the null-result, wh ile iRegulo fc Roca Cortes! (|2005l ) and 
iMarchenkol (|2008l ) cautiously claimed a detection of os- 
cillati ons based on a re-an alysis of MOST data. Mean- 
while, iBruntt et al.l (|2005l) reported a detection of power 
excess with amplitudes of about twice the solar value (~ 
8 ppm) based on continuous space-based photometry by 
the WIRE satellite, consistent with the upper limits set 
by the MOST results. 

The clear identification of individual oscillation modes 
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in Procyon was finally achieved with a large ground- 
based radial velocity campaign that was carried out i n 
January 2007 (jArentoft et a l. 2008; B edding et aHl201~0l ). 
Simultaneously, a third set of MOST observations, longer 
and with higher precision than the pre vious runs, was 
obtai ned in January and February 2007 (jGuenther et al.1 
2008). Here we present the first direct comparison of 
these datasets. 

2. SUMMARY OF OBSERVING CAMPAIGNS 

The following sections present a brief summary of the 
observations and main results of the two campaigns on 
which the comparison in this paper is based. For a brief 
introduction to basic characteristics of solar-like oscilla- 
tion s, relevant to P r ocyon , we refer the reader to Section 
2 in lBedding eTail (|2010t ). 

2.1. MOST photometry 

The MOST (Microvariability and Oscillations of 
STars) space telescope, launched in 2003, is the first 
satellite dedicated to asteroseismic observations from 
space. MOST houses a 15 cm telescope, with obser- 
vations performed through a custom broad-band filter 
(350-700 nm) . It is positioned in a sun-synchronous low- 
Earth orbit, enabling it to continuously monitor stars 
for up to three months. For a detailed desc ription of the 
i nstrument, we re fe r tolWalker et all (|2003l ). 

iGuenther et al.l (|2008l ) presented the results of the 
MOST 2007 campaign. The data are the most pre- 
cise MOST Procyon photometry to date, outperform- 
ing the previous runs in both length (38.5 d) and time- 
series point-to-point scatter (140 ppm), and consequently 
also in high-frequency noise (0.9 ppm). While the power 
spectrum showed an excess in the expected oscillation 
frequency range and an autocorrelation of the spectrum 
yielded strong evidence for the expected characteristic 
large frequency separation (Af) of ~ 55/iHz, no reliable 
individual mode frequencies could be extracted. Us- 
ing various common a nd new frequency extraction tools, 
IGuenther et al.l (j2008fl found that no consistent regularly 
spaced frequencies could be identified unless it was as- 
sumed that the mode lif etimes are ^2 days, i.e. slightly 
less than solar (see, e.g.. iChaplin et aLll2009h . The large 
scatter of the extracted frequencies around the predicted 
regularity (±5 /xHz) was found to be consistent with the 
scatter of the frequency identifications from previous ra- 
dial vel ocity campaigns. Using new numerical convection 
models IGuenther et al.1 (|2008l ) argued that, contrary to 
the Sun, the granulation timescale in Procyon is simi- 
lar to the timescale of the p-mode oscillations. This was 
identified as a possible reason for the short mode lifetimes 
in Procyon and the consequent difficulty of extracting 
consistent frequencies from the MOST 2007 data. 

2.2. Ground-based radial velocity campaign 

The 2007 ground-based radial velocity (RV) c ampaign 
was described in detail bv lArentoft et all (|2008l ) and in- 
cluded 11 telescopes, with apertures ranging from 0.6 to 
3.9 m, at eight observatories. Covering a total length of 
25 days, with a duty cycle above 90% for the central 10 
days of observations, it is the most precise and complete 
radial-velocity campaign dedicated to asteroseismology 
to date. 



lArentoft et al.1 (|2008l ) reported the detection of slow 
variations in the velocity timeseries with an apparent pe- 
riod of ~10d, which they attributed to rotational mod- 
ulation of active regions on the surface of the star. A 
possible rotation period of Procyon of ^10 d or twice 
that value was suggested, with the latter scenario be- 
ing more likely if it is assumed that the rotation axis 
is al igned with the known inclination of the binary or- 
bit. lArentoft et al.l (|2008l ) also provided an estimate of 
the mean mode amplitude between 650 — 1150/iHz of 
38.1±1.3 cms -1 , consistent with previously reported de- 
tections and upper limits. 

A detailed asteroseismic a nalys is of the RV data was 
presented bv lBedding et al.l (|2010ft . The continuous cov- 
erage and high S/N allowed the first measurements of 
mode frequencies in Procyon. These authors also mea- 
sured large and small frequency separations as a func- 
tion of frequency and identified a possible mixed mode 
at 446 /iHz. They used the variation of peak amplitudes 
caused by the stochastic nature of the oscillations to es- 
timate the mode lifetimes in Procyon to be 1.3±0.5d. 
As a result of the large linewidths of the modes, how- 
ever, a major difficulty in the analysis was the identi- 
fic ation of modes with even and odd spherical degree 
I. iBedding et al.l (|2010l ) gave several arguments for the 
most probable mode identification but left open the pos- 
sibility that the identification could be reversed. 

3. SLOW VARIATIONS 

The Procyon MOST photometry provides the oppor- 
tunity to further investigate the nature of the slow varia- 
tions detected in radial velocity. The stability and conti- 
nuity of MOST photometry has previously b een success- 
fully used to study activ ity in several stars (jCroll et al.l 
IliMlWalker et al.|l200l. 

The analysis bv IGuenther et al.l ((2008) was based on 
a high-pass filtered MOST light curve to focus on the 
detection of solar-like oscillations. Here, we use the raw 
light curv e, as produced by t he reduction pipeline de- 
scribed in iReegen et al.1 ()2006[ ) , as the start i ng poin t of 
our analysis. As discussed by iHuber et al.1 (|2009b[ ) for 
another solar-like star observed by MOST (85 Peg) the 
reduced photometry sometimes shows long-periodic in- 
strumental variability that can be identified and cor- 
rected by decorrelating satellite telemetry data, such as 
board and preamplifier temperatures, against the ob- 
served target intensities. We applied the same technique 
to the 2004, 2005 and 2007 photometry of Procyon and 
the resulting detrended light curves are shown in Figure 

The intensity and velocity curves in Figure[T]show sim- 
ilar variability in all datasets. The most prominent peri- 
odicity in the velocities, with a period of ~10 days, is not 
readily apparent in the 2007 MOST photometry, which 
shows the strongest signal as measured from the ampli- 
tude spectrum with a period of ~6 days. How ever, as 
is wel l known for the Sun and demonstrated by iClarkej 
(2003) for models of spotted stars, the relationship be- 
tween simultaneous velocity and intensity observations 
of active regions is not simple and strongly depends on 
parameters such as rotational velocity, inclination and 
spot size. Detailed modeling of the variations using the 
overlapping parts of the dataset (~ 10 days) is beyond 
the scope of this paper, but we make some qualitative 



Oscillations and activity in Procyon 



3 



a 1000 

ft 

Q. 

~-~ 500 

a 

£ -500 

« -1000 



^^^^^^^ 



MOST 3004 



i oooo 



1480 1490 1500 

HJD- 345 1545 (d) 




1870 I860 

HJD -345 1545 (d) 



1000 
500 


-500 
-1000 

3000 
1000 



-1000 
-3000 





Mi 




; MOST 3007 


; *B' ."1 


♦ *? 






! • s 

1 ** tkii •> 






EV 3007 1 




Pj5f 







3560 3570 2580 3590 3600 

HJD -345 1545 (d) 

Fig. 1. — Detrended MOST light curves of Procyon obtained in 
2004, 2005 and 2007 (every 50th datapoint shown) and 2007 radial 
velocity curve (every 20th datapoint shown) of Procyon. The thick 
black lines show the result of smoothing with a boxcar of width ~ 
2 days. 



statements on the variability based on the observed low 
frequency power levels. 

To compare both datasets independent of length and 
sampling, we converted the power spectra to power den- 
sity by multiplying with the effective length of each 
dataset (calculated as the inverse of the area under the 
spectral window in power). Figure [2] shows the low- 
frequency power density level measured in Procyon in 
intensity (upper panel) and velocity (lower panel) using 
the full length of both datasets obtained in 2007. To 
measure the low frequency power density level, we fitted 
power laws with a fixed slope of 2 in the frequency in- 
terval 5-40 /iHz (red dashed lines). For comparison, the 
blue dashed-dotted lines show the average power density 
level of the Sun during solar maximum derived using the 
same method with 30 d subsets of data obtained in in- 
tensity by VIRGO (iFrohlich et al]|1997ft and in v elocity 
by GOLF (lUlrich et all 120001 IGarci'a et alJl2005h . both 
of which are instruments onboard the SOHO spacecraft. 
Note that we have corrected the photometric power den- 
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Fig. 2. — (a) Photometric power density spectrum of Procyon. 
The red dashed line shows a power law with a fixed slope of 2 
fitted in the frequency interval 5—40 /iHz. The blue dashed-dotted 
line shows the average power density level of the Sun during solar 
maximum, determined using the same method with 30 d subsets of 
SOHO data, (b) Same as panel (a) but for radial velocities. 



sity for the spe ctral response of the MOST filter following 
the method of iMichel et~aTl (|2009t ). yielding R g = 4.23 
using an effective temperature of T c g — 6500 K for Pro- 
cyon, compared to R g = 5.02 for solar VIRGO observa- 
tions in the gr een channel. 

As noted bv lArentoft et all (|2008[ ). the Procyon veloc- 
ity power spectrum shows a similar 1/v 2 dependence as 
observed for the Sun but at higher power density lev- 
els, and the same is observed in the photometry. Is the 
observed excess of photometric power density, measured 
relative to the Sun, consistent with the velocity obser- 
vations? To investigate this, we used SOHO data span- 
ning from 1996 to 2004 and measured the low frequency 
power density levels in independent 30 d subsets as de- 
scribed above throughout the solar activity cycle. The 
ratios of the observed levels between Procyon and the 
Sun for each subset are shown in Figure EJa). 

We observe that the power density ratio between Pro- 
cyon and the S un is substantially hig her in velocity than 
in photometry. lArentoft et al.l ()2008l ) argued that the ve- 
locity power density at frequency v is expected to scale 
as follows: 



, da 

PD {v RV = — 
1 a 



v sin i 
T 



(1) 



where da /a is the fractional area covered by active re- 
gions, vsini the projected rotational velocity, and T the 
typical lifetime of active regions on the stellar surface. 

Similarly, the amplitude measured in photometry will 
be proportional to da/a and the luminosity variation 
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SL/L caused by the flux contrast between the unspotte d 
and spotted areas of the star (see, e.g., lDorrenlfi987l) . 
We can therefore rewrite Equation (1) for the case pho- 
tometric power densities as: 



PD(i/) 



Phot 



(2) 



We note that Equation (2) is only intended to give an 
approximate estimation of the photometric power den- 
sity due to stellar activity. C ompared to detailed photo- 
metr i c spot models (see , e.g. lDorreniri987t lLanza et al.l 
120031: iMosser et alj|2009h . Equation (2) for example does 
not include a n explicit depe n dence on stellar inclination. 
As shown by IMosser et al.l fj2009f ) . this corresponds to 
neglecting any information about the latitude A of ac- 
tive regions (i.e., we assume in Equation (2) that it is 
equally likely to observe spots near the equator or near 
the poles). While this is certainly not the case for the 
Sun, this information is also neglected in Equation (1) 
(e.g., an active region at A = 90° on a star with i = 90° 
will cause no velocity variation). Since we are here only 
interested in comparing the ratio of velocity to photom- 
etry variations, any explicit dependence on the latitude 
of active regions will therefore cancel out. 

Equations (1) and (2) imply that the difference be- 
tween the velocity and photometry power density ratio 
of two stars depends only on usini and SL/L. We as- 
sume that the ratio in SL/L (and hence the flux contrast) 
between Procyon and the Sun is negligible compared 
to the ratio in usini, which is supported by detailed 
spot modeling of stars hotter than the Sun (jLanza et al.1 
120091 120111 ). Using vsini = 2.0 km s" 1 for the Sun and 
v sin i = 3.2 km s -1 for Procyon (|Allende Prieto et al.l 
2002), this therefore implies a difference in the veloc- 
ity and photometry power density ratios by a factor of 
~ 2.56. The red dashed-dotted line in Figure [3ja) shows 
the photometric power density ratio multiplied by this 
value. We observe that the power density ratios are 
now in better agreement, both during solar minimum 
(in ^1996) and solar maximum (in ^2002). This result 
implies that the observed variations of Procyon in both 
datasets are qualitatively in agreement with being the 
signature of stellar activity. 

The scatter of individual datapoints in Figure [3ja) 
is large, presumably due to the complex non-sinusoidal 
variations of stellar activity causing large variations in 
low frequency power density fits. We have therefore re- 
peated t he above exercise u sing the activity range as de- 
fined in Bas ri et al.l (|2010| ). which measures the maxi- 
mum absolute deviation of the time series with respect 
to its mean. The results of this are shown in Figure 
[2Kb). In this case, the scaling factor is expected to be 

V2.56. Again, the scaled photometric ratios are in rea- 
sonable agreement with the velocity ratios, confirming 
the results derived using power densities. 



4. SOLAR-LIKE OSCILLATIONS 

4.1. Comparison of power spectra 

Figure 0] compares the power spectra of both full 2007 
datasets in the frequency range where p modes have been 
detected. We used the sidelobe-optimized power spec- 
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Fig. 3. — (a) Ratio of low-frequency power density for Procyon 
and the Sun as a function of solar activity cycle in photometry 
(triangles, dashed line) and velocity (diamonds, solid line). Open 
symbols show individual measurements using 30-day subsets and 
lines are values which were smoothed twice using a boxcar with a 
width of ~5 months. The red dashed-dotted line shows the photo- 
metric ratio scaled to account for the different projected rotational 
velocities of Procyon and the S un, (b) Same as Pa nel (a) but using 
the activity range as defined in Basri et al. (2011]). 



tram for the velocities, as described by Bedding et al. 
(12010ft. and the high-p ass filtered MOST data used by 
IGuenther et al.l (|2008f ). Note that the funnels of low 
power around the orbital harmonics of MOST (marked 
as dashed lines) are due to the high-pass filter remov- 
ing power leaking from low frequencies, which in turn is 
caused by periodic outlier rejections in the data reduc- 
tion pipeline during high-straylight phases. To smooth 
over the effects caused by the stochastic nature of the os- 
cillations, we convolved both spectra with a Lorentzian 
profile with a width of 2.5/dTz, corresponding to a 
mode lifetime of 1.5 d. We also show with dotted lines 
the values for the odd and even ridge centroids deter- 
mined by Beddin g et al.1 ((2010) for the velocity data 
(see their Figure 9). Note that although our adopted 
value for the mode lifetime is larger than measured for 
other F-stars su ch as HD 49933 (jGruberbauer et alj2009t 
iBenomar et al.|[2009f ). the exact choice of this value has 
no influence on the results presented below. 

The comparison shows clearly that most of the peaks 
in the MOST spectrum coincide with pulsation frequen- 
cies identified in the velocity power spectrum. The 
agreement improves considerably towards high frequen- 
cies (> 800/xHz), where the granulation background in 
photometry becomes lower, as can be seen by the steady 
decrease of power towards high frequencies (see also Sec- 
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Fig. 4. — Power spectra of Procyon from MOST photometry (upper panel) and the radial velocity campaign (lower panel), smoothed 
with a Lorcntzian with a width of 2.5/iHz (corresponding to a mode lifetime of 1.5 d). Dotted lines mark the odd and even ridge centroids 
identified in the velocity data. Dashed lines show the harmonics of the MOST orbital frequency. 



tion 14. 3p . We also see that, unluckily, several of the in- 
trinsic pulsation frequencies of Procyon near maximum 
power coincide almost exactly with harmonics of the or- 
bital frequency of the MOST satellite. 

To investigate the agreement between the two spec- 
tra in more detail, we calculated a cross-correlation of 
the two power spectra in the central region of maximum 
power from 650-1150 /iHz. The result is shown in Fig- 
ure [5l As expected, we see clear maxima at zero offset 
and at half the large frequency separation. Note that the 
shift of the highest peak from zero offset is only 0.9 /iHz, 
which is much smaller than the lifetime of the modes 
and therefore insignificant. To test the significance of 
the peak height, we correlated the RV power spectrum 
with MOST power spectra calculated from 2000 white 
noise time series with the same sampling and scatter as 
the original dataset. The resulting distribution at zero 
offset showed a mean of zero and a standard deviation 
of 0.09, which is indicated as dashed and dotted lines 
in Figure [5] At a level of > 4er, these results confirm 
that the peaks observed in the MOST power spectrum 
very likely correspond to the oscillations observed in the 
radial velocity data. 

4.2. Echelle diagrams and folded power spectra 

A widely used method to analyze the regular frequency 
pattern characterizing solar-like oscillations is to stack 
the power spectrum (or extracted frequencies) in slices 
of the large freq uency separation At/, forming a so-called 
echelle diagram (jGrec et al.lll983l ) . Note that throughout 
the paper we use Av — 56 /iHz, which corresponds to the 
large separation at maximum power as identified in the 
RV data (see Figure 11(a) in lBedding et^[2010t) . The 
left and middle panel of Figure [6] show echelle diagrams 
of both power spectra smoothed to the same frequency 
resolution. Note that the MOST power spectrum has 
been corrected for the background contribution due to 
granulation and activity (see Section |4~5|) . The echelle 
diagrams clearly show two ridges corresponding to modes 
of odd and even degree in both datasets. The similarity 
of the curvature of both ridges (and hence the variation 
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Fig. 5. — Cross-correlation of the velocity and photometry power 
spectra in the frequency range 650-1150 /iHz. The dashed and dot- 
ted lines mark the mean and ±lcr noise level derived from cross- 
correlating the RV power spectrum with 2000 MOST power spectra 
computed from white noise time series with the same sampling and 
scatter as the original data. 



of Av with frequency) in both individual datasets reaf- 
firms our conclusion that the peaks seen in Figure U are 
intrinsic p modes. The echelle diagrams also show that 
the higher low-frequency noise in the photometry makes 
it harder to detect p modes below ~ 800/iHz than in the 
velocity data. 

In addition to the MOST and RV datasets, we also 
analyzed a power spectrum constructed by multipliying 
both individual power spectra. This corresponds to a 
power spectrum of the convolution of the MOST and 
RV timeseries, combining the advantage of the higher 
frequency resolution of the MOST data with the higher 
S /N of the RV data. The resulting echelle diagram dis- 
played in the right panel of Figure |5] clearly shows more 
well-defined ridges than the corresponding diagrams us- 
ing the individual datasets. 

Having confirmed that we have detected the ridges in 
both datasets, an obvious step in the analysis is to at- 
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Fig. 6. — Echelle diagrams of MOST (left panel) and radial velocity (middle panel) power spectra, as well as the combination of both 
data sets (right panel). Darker greytones correspond to higher power. The asterisks in the left and right panel mark the harmonics of the 
MOST orbital frequency. 



tempt to confirm or discard th e ridge identification pre- 
sented in lBedding et all (|2010l ) using the MOST data. It 
is well known that observations in intensity are less sen- 
sitive to pulsation modes of higher spherical degree than 
velocity measurements. We therefore expect ridges ob- 
served in photometry to be shifted to higher frequencies 
in the echelle diagram than in velocity, with the amount 
of shifting depending on the ridge identification. Cal- 
culations based on th e theoretical response functions by 
IKieldsen et al.l (2008) and the frequency values presented 
in IBedding et al] (]2010f ). however, showed that the ex- 
pected shift is only of the order of 1 /LtHz, too small to be 
detected with the current uncertainties. 

Another possibility to search for ridge asymmetries is 
to increase the S/N by collapsing the echelle diagram 
over several orders. To do so, we calculated the ridge 
centroid s for the MOST data in the same manner as 
done by IBedding e t al.l (|201Q[ ) for the radial velocities. 
Using the ridge centroids calculated for the MOST data, 
we straightened each order by removing the curvature 
seen in the echelle diagram before collapsing the echelle 
diagram in the frequency range where the centroids were 
reliably determined (700-1300 /iHz). Figure [7] shows a 
comparison with the analogous pro cedure for the veloc- 
ity power spectrum (see Figure 10 in lBedding eFaLll^oToT) 
as well as for the product of the MOST and radial veloc- 
ity power spectrum. Note that the ridge centroids have 
been calculated for each of the three power spectra in- 
dividually since, as explained above, it is expected that 
the ridge positions are slightly different for each dataset. 
While the two ridges are clearly detected in the MOST 
photometry (top panels), the S/N is too low to make 
any firm conclusion about the possible presence of sepa- 
rated ridges (i.e. I = and 1 = 2), as is the case for the 
RV data (middle panels). We note, however, that the 
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Fig. 7. — Collapsed echelle diagrams in the region 700 — 1300 fiHz 
after correcting for curvature using ridge centroids. The top panels 
show the diagrams calculated using the MOST power spectrum, 
the middle panels using the radial velocity power spectrum, and 
the bottom panels using the product of both power spectra. Panels 
are separated into the left-hand side ridge (panels a, c and e) and 
right-hand side ridge (panels b, d and e) seen in Figure |S] The 
dotted lines in the middle panels show the radial velocity power 
spectrum colla psed over the full osc illation range (corresponding 
to Figure 10 in lBedding et al.l ||20T0T )') 



collapsed power spectrum of the combined data (bottom 
panels) might show some evidence for separated I = 
and I = 2 components in the right-hand-side ridge. This 
would im ply Scenario A, which is opposite to that pre- 
ferred by Bedding et al] (|2010[) but in agreement with 
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the r esults of Bayesian model comparisons using the RV 
data (jBedding et alJ2010tlHandberg fc Campantdl2011[ ). 
A more detailed analysis including the extraction of indi- 
vidual frequencies from the combined data and a compar- 
ison with pulsation models will be presented in a forth- 
coming paper (T. Kallinger et al., in preparation). 

4.3. Oscillation amplitudes 
4.3.1. Amplitude Ratios 

The simultaneous observing campaigns allow us to 
measure the ratio of oscillation amplitudes in photome- 
try and velocity and therefore test t heoretical values and 
the sc aling relations introduced by iKieldsen fc Bedding] 
(1995). The following influences have to be considered 
when measuring amplitudes of solar-like oscillations: 

(i) Systematic variations due to the stellar cycle effects 

(ii) Variations due to stochastic excitation and damp- 
ing of the oscillation signal 

(iii) Measurement errors due to background subtraction 

Since both datasets have been obtained within a times- 
pan of less than 60 days, effects arising from (i) can be 
safely ignored here. To measure amplitudes in a way 
that is largely insensitive to point (ii), we convolved the 
power spectrum with a Gaussian with FWHM = 4Av 
and scaled the signal to the contribu tion of radial modes 
in each order (|Kieldsen et al.l 120081 ) . Note that for the 
factor c, which measures the effective number of modes 
per order, we have i n terpo lated the values listed in Table 
1 of IKieldsen et al] (|2008l ) to the central wavelength of 
the MOST filter (525 nm), yielding c = 3.14. 

A crucial precondition to estimate the amplitude, par- 
ticularly at low S/N, is to properly correct for the back- 
ground contribution arising from stellar granulation and 
activity. In order to reliably estimate the uncertainty 
of the background parameters, we have used a combina- 
tion of two published m ethods. An initial le ast-squares 
fit using the method of iHuber et alj ((200 9a) was used 
as a starting point for a more careful fitting procedure 
using a Bayesian Markov -Chain Monte-Carlo (MC MC) 
algorithm, as describe d in lGruberbauer et al.l (|2009| ) and 
iKallinger etTaD (120101) . The fitted model was adopted 
fromlKarofl ( 2008) and has the form: 

where P n is the white noise component, k is the number 
of power laws used and a and t are the rms intensity and 
timescale of granulation, respectively. Note that in our 
application for Procyon, k = 2. 

The determination of the resulting background and 
corresponding amplitude was done in two steps: 

1) We excluded the region of the power spectrum 
which contains oscillation signal and ran the 
Bayesian MCMC algorithm only on the remain- 
ing power spectrum. The excluded region was de- 
termined by visual inspection of heavily smoothed 
power spectra (see Figure [8]), and we verified 
through several trial MCMC runs that changing 
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Fig. 8. — Power density spectra of Procyon smoothed with a 
1 /^Hz boxcar (light grey) and 110 fiHz boxcar (dark grey) observed 
in photometry (upper panel) and velocity (lower panel) with fitted 
background models (solid black lines). Vertical dotted lines mark 
the frequency intervals occupied by stellar oscillations which were 
excluded from the fit. The insets show the distributions of the mean 
amplitude per radial mode in the frequency interval 650-1150 (j,Hz 
derived from the MCMC simulations of the background fit. Dashed 
lines mark the median and dotted lines the 1-cr confidence limits. 

this region within reasonable limits does not sig- 
nificantly influence the results. We used uniform 
priors for the granulation timescales and Jeffreys 
priors for the amplitudes of the background compo- 
nents. The most probable background model was 
then determined as the median of the marginalized 
posterior distributions for each background param- 
eter. 

2) Since the amplitude is not implicitly defined as a 
model parameter that is fitted to the data, we de- 
termined the pulsation amplitude as the mean level 
of the smoothed power spectrum in frequency range 
650-1150 /iHz after subtracting the background. 
This step was done for each MCMC iteration, and 
the final amplitude was then evaluated as the me- 
dian and l-cr confidence limits of the resulting am- 
plitude distribution. 

Figure [5] shows the background fits resulting from 10 
independent MCMC chains, each with 10 iterations, 
for the full MOST and radial velocity datasets. The 
distributions for the mean amplitude over the range of 
650-1150 /iHz (after subtracting the background signal) 
yields Ai =0 , p hot = 9.1Io;4Ppm for the MOST data and 
^2=0. rv = 40.0^04 cm -1 s for the velocity data. Hence, 
we arrive at a mean amplitude ratio for the full datasets 
of A=n. pr ,nt.A4i =o,Ry = 0.23 ± .01 pp m cm -1 s. 

As shown by IKieldsen et al.l (|2008| ) for the Sun, short 
mode lifetimes can cause considerable variations even 
after heavily smoothing the oscillation envelope. To 
test these effects, we employed a similar approach as 
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lArentoft et all (|2008[ ) and subdivided the MOST data 
into five independent subsets of equal length (~ 7.5 d). 
Note that this subset length is considera bly longer than 
the 1 individual 2 d subsets used by lArentoft et al.1 
( 2008) for the radial velocity data since the signal in pho- 
tometry is lower, requiring longer sets to achieve a suffi- 
ciently high signal for an amplitude determination. The 
resulting smoothed curves, which were corrected for a 
background ht calculated using 5 x 10 5 MCMC iterations 
for each subset, together with their mean value are shown 
in Figure Eta). Using the central region of maximum 
power between 650-1150 /^Hz, we derived a mean ampli- 
tude of A/ = o,phot = 9.1 ± 0.5 ppm, in very good agree- 
ment with the value derived using the full dataset above. 
In order to compare this photometric value to a velocity 
amplitude derived using the same method, we combine it 
wi th the estimate for the radial velocity amplitude given 
bv lArentoft et all (|2008l ) as ^ =0 ,rv = 38.1 ± 1.3 cm- 1 s. 
This yields an amplitude ratio, measured using subsets, 
of A; =0 ,photA4i=o,RV = 0.24 ± 0.02 ppm cm -1 s, again 
in good agreement with the value derived from the full 
datasets. 

A few points need to be considered when evaluat- 
ing the quoted values and uncertainties. Firstly, the 
fact that the uncertainty on Ai = o t nv derived from the 
full dat aset is considerably l ower than the value de- 
rived by lArentoft et ail (|2008| ) shows that for high S/N 
data the influence of stochastic excitation (see point (ii) 
above) dominates the amplitude uncertainty over the 
uncertainty arising from the background determination 
(see point (hi) above). For the MOST data on the 
other hand, both uncertainties are roughly the same. 
To ensure a conservative approach, we therefore opted 
to use the amplitude ratio determined using subsets, 
^4;=o,photM;=o,RV = 0.24 ± 0.02 ppm cm -1 s, as our fi- 
nal value of the mean amplitude ratio. 

Secondly, the amplitude ratios do not include any 
uncertainties arising from the factors c used to nor- 
malize the amplitude per radial mode. First observa- 
tio nal constraints o n phot ometric mode visibility ratios 
by iDeheuvels et al.l (|2010| ) suggest ~ 2a differences of 
up to 25% compared to theoretical responses used by 
iKieldsen et al.l (|2008l ). These differences could be caused 
by uncertainties in the limb-darkening laws used to cal- 
culate theoretical response functions, but also intrin- 
sic differences in amplitudes of different degrees. We 
will assume that the latter cancel out in the photome- 
try to velocity ratio, and therefore only concentrate on 
the mode visibilities. We repeated the calculations of 
c by numerically integrating the spatial response func- 
tions for a quadrati c limb-darkening law for the Sun 
(jBedding et al.l Il996l ) assuming a conservative absolute 
uncertainty of 0.05 (corresponding to a relative uncer- 
tainty of ^10- 20%) for each limb-darkening coefficient 
(|Howarthll20lot) . The resulting c factors after 1000 inte- 
grations yield an uncertainty of 4% in velocity and 6% 
in photometry, which translates into uncertainties of 2% 
and 3% in the normalized amplitudes, and hence an un- 
certainty of 4% in the amplitude ratio. This test shows 
that the uncertainty on c can be substantial for estimat- 
ing amplitude ratios, in particular for observations with 
high S/N and long mode life times as found in red giant 
stars (jDe Ridder et al.l [20091 ) . In our case, however, the 
amplitude ratio is dominated by the ~8% uncertainty 



arising from the background fits and finite mode life- 
times. 

4.3.2. Comparison with Theoretical Results 

How do our estimates for the amplitude ratio com- 
pare with scaling rel a tions? Rearranging equation (5) in 
IKieldsen fc Beddind (fl995l) yields 

(6L/L) x /ppm = 20.1 x HT 3 

^osc/cms- 1 (A/550nm)(T cff /5777K)'' ' 1 ' 

with r = 1.5 if the oscillations are adiabatic, and a best- 
fitting value o f r = 2 for observed amplit udes in classical 
pulsators (see IKieldsen fc Bcddingl ll995() . 

Using r = 2, T c g — 6500 K and a wavelength A = 
525 nm, the expected amplitude ratio for radial modes is 
0.17ppmcm~ 1 s. As shown by the dotted line in Figure 
Eld), this is considerably lower than the observed ratio. 
Using r — 1.5 the agreement is only slightly better, with 
an expected ratio of 0.18 ppm cm" 1 s. Since the radial 
velocity amplitude is in agreement with t he value from 
scaling relations (see lArentoft et al.l l2008f) , this implies 
that the MOST amplitude is higher than expected. 

Going one step further, the smoothed amplitude curves 
allow us to analyze the amplitude ratio as a function of 
frequency. Figure E^d) shows the ratio of the amplitude 
curves derived from the full datasets for MOST and the 
RV data, which are shown separately in Figure Eljb) and 
Etc). Note that we have scaled the uncertainty of the 
RV amplitude curve to match the relative uncertainty 
derived using subsets, as described in the previous sec- 
tion. We also repeated the calculation using only HARPS 
data from the RV dataset, and restricted the MOST data 
to the same timespan (~6.5d). The result was almost 
identical to the result based on the full dataset but with 
larger uncertainty. 

Figure Ete) compares the amplitude ratio as a func- 
tion o f frequency with theoretical predictions by iHoudekl 
(2010) for a model of Procyon for different model at- 
mospheres and different heights above the photosphere. 
Although the exact shape of the variation as a function 
of frequency is not well recovered, we note that the aver- 
age amplitude ratio at maximum power (~ 1000/iHz) is 
in better agreement with models than the estimate based 
on scaling relations (dotted line). A more detailed com- 
parison will have to await the collection of higher S/N 
data from which amplitudes of individual mode frequen- 
cies can be reliably extracted. 

It is interesting to note that, as can be seen from 
Figures EJb) and Efc), the maximum of the os- 
cillation envelope in photometry seems to be shifted to 
slightly higher frequencies than in velocity. Defining u max 
as the frequency corresponding to the maximum of the 
smoothed oscillation envelope, the MCMC analysis of the 
full datasets yields f m ax,phot = 10141^/iHz, compared to 
^max,RV = 9231^/iHz. This shift translates into an in- 
crease of the amplitude ratio as a function of frequency, 
which tentatively can be identified in Figure EJd) . 

5. CONCLUSIONS 

We have compared simultaneous space-based MOST 
photometry and ground-based radial velocity data of the 
F5 star Procyon. Our main findings can be summarized 
as follows: 
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• The peaks observed in the MOST 2007 power spec- 
trum match the oscillation frequencies detected in 
the radial velocity campaign, and the echelle dia- 
grams shows similar structure and curvature. We 
have attempted to confir m the mode r idge i den- 
tification presented by iBedding et al.l ()2010D by 
collapsing power spectra corrected for curvature. 
While the MOST data alone do not provide con- 
clusive results, the collapsed echelle diagram of the 
combined datasets shows some evidence that sce- 
nario A in lBedding et al.l (|2010l ) is the correct mode 
identification, contradicting the conclusions of that 
paper. Further work based on the extraction of in- 
dividual frequencies from the combined data and 
a comparison with stellar models will be necessary 
to confirm this result. 



• We measured the mean luminosity amplitude per 
radial mode in Procyon in the frequency range 
650-1150 /^Hz to be Ai =0tP ^ ot = 9.1 ± 0.5 ppm, 
in agreem ent with the v alue o f 8.5 ± 2 ppm pub- 
lished by iBruntt et all (|2005| ). Combining this 
with the mean veloc ity amplitude measured by 
lArentoft et al.l (|2008t ) gives an amplitude ratio 
of A i= o !P hotMi=o,RV = 0.24 ± 0.02 ppm cur 1 s. 
This is considerably higher than the value of 
0.17 ppm cm -1 s expected from scaling from the 
Sun, but is in bett er agreement w ith theoretical 
values predicted by|Hpudck ( 20101) . We also ana- 
lyzed the amplitude ratio as a function of frequency 
and found that the maximum of the oscillation en- 
velope appears to be shifted to higher frequencies 
in photometry than in velocity. 



The results presented here illustrate the potential of 
combining simultaneous luminosity and velocity mea- 
surements to study pulsations in stars. Future op- 
portunities may arise from combining measurements 
from the ground-base d radial-velocity network SONG 
(|Grundahl et al.l 120081) with space-bas ed photometry by 
MOST and BRITE (| Weiss et all [20081 ) . which will mark 
an important step in studying stellar structure and evo- 
lution in bright stars with well-known fundamental pa- 
rameters. 



Fig. 9. — (a) Smoothed amplitude curves of five independent 
7.5 d subsets of MOST photometry (thin lines) together with their 
mean, (b) Smoothed amplitude curve using the full MOST datasct. 
Dashed lines mark the 1-<t uncertainties, (c) Same as panel (b) but 
for the full RV timeseries. (d) Ratio of the curves shown in panels 
(b) and (c). The horizontal dotted lin e shows the expected ampli- 
tude ratio using the scaling relation of Kieldsen & Bedding ( 1995). 
(e) Same as panel (d), but compared to theoretical amplitude ra- 
tios using an Eddington atmosphere (red lines) and a scaled VAL-C 
atmosphere (blue lines), for heights of 300 km (dashcd-dotted lines) 
and 600 km (d ashed-triple-dotted lines) above the photosphere (see 
Houdek 2010, for details). 

• The MOST light curves of 2004, 2005 and 2007 
show slow variations similar to those observed in 
velocity. A comparison of the variability level in 
photometry and velocity to the Sun confirmed that 
these variations are compatible with stellar activity 
on Procyon. 
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